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ABSTRACT 

We perform direct numerical simulations of the tidal encounter of a rotating 

s^ , planet on a highly eccentric or parabolic orbit about a central star formulated as an 

^ '. 

r^ , initial value problem. This approach enables us to extend previous work of Ivanov 

0\ 
CO ■ & Papaloizou to consider planet models with solid cores and to avoid making an 

(N 

o 
o 



anelastic approximation. 
^D I We obtain a power spectrum of the tidal response of coreless models which en- 



ables global inertial modes to be identified. Their frequencies are found to be in good 
agreement with those obtained using either a WKBJ approach or the anelastic spec- 

^% ' tral approach adopted in previous work for small planet rotation rates. We also find 

^ '■ 

5^ , that the dependence of the normal mode frequencies on the planet angular velocity 

in case of higher rotation rates can for the most part be understood by applying first 

order perturbation theory to the anelastic modes. 

We calculate the energy and angular momentum exchanged as a result of the 
tidal encounter and for coreless models again find good agreement with results ob- 
tained using either the anelastic spectral method. 

Models with a solid core showed evidence of the emission of shear layers at 
critical latitudes and possibly wave attractors after the encounter but the total energy 
exchanged during the encounter did not differ dramatically from the coreless case as 
long as the ratio of the core radius to the total radius was less than 50%, there being 
hardly any difference at all when this ratio was less than 25% of the total radius. We 
give a physical and mathematical interpretation of this result. 

Finally we are able to validate the use of the anelastic approximation for both the 
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work presented here and our previous work which led to estimates of circularisation 
rates for planets in highly eccentric orbits. 

Key words: hydrodynamics; stars: oscillations, binaries, rotation; planetary sys- 
tems: formation 

1 INTRODUCTION 

A significant number of extrasolar planets are found in close circular orbits around their central star 
while others are found in highly eccentric orbits. This is an indication that planet-planet scattering 
may have been important after formation and generated highly eccentric orbits that subsequently 
undergo circularisation with attendant tidal heating ( see Papaloizou & Terquem 2006 for a review 
and references therein). In order to investigate such a scenario the tidal interaction of a planet on 
a highly eccentric orbit with the central star has to be evaluated. This was considered in the non 
rotating case by Ivanov & Papaloizou (2004) and in the rotating coreless case using an anelastic 
approximation by Papaloizou & Ivanov (2005) and Ivanov & Papaloizou (2007) (hereafter referred 
to as PI and IP, respectively). The excitation of pulsations or normal modes of the planet is a key 
aspect of the tidal problem. 

In this paper we extend our previous work by relaxing the anelastic approximation and con- 
sidering planet models with solid cores. We investigate pulsations of a uniformly rotating fully 
convective or barotropic object by performing numerical simulations of tidal encounters consid- 
ered as initial value problems. The formalism we adopt does not require an anelastic approximation 
to be made and so it enables both the response of the low frequency inertial modes of oscillation 
as well as the higher frequency / and p modes to be taken into account. As in our preceding stud- 
ies we consider the pulsations to be excited in an object referred to hereafter as a 'planet' as a 
consequence of moving on a parabolic orbit around a central object. 

We calculate the energy and angular momentum exchanged after pericentre passage for a va- 
riety of pericentre distances and angular velocities of rotation of the planet. We also study the 
spectrum of normal modes that is excited and compare their properties with the analytic theory 
developed in our accompanying paper Ivanov & Papaloizou (2009) hereafter IPN and also results 
obtained using the anelastic spectral approach described in Ivanov & Papaloizou (2007) (here- 
after referred to as IP) . This comparison finds good agreement between the different approaches, 
validating the use of the anelastic approximation for describing the inertial mode response. 

We also consider the effects of introducing a solid core into the planet. Although we find 
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evidence for the development of singular phenomena such as the emission of shear layers at 
critical latitudes and the formation of wave attractors (see Ogilvie & Lin 2004, Ogilvie 2009 and 
Rieutord & Valdettaro 2010 and references therein) and after the encounter, the energy exchanged 
is in general not changed by a very large amount and for a modest sized core hardly at all. We give 
both a physical and mathematical explanation for this result in appendix B. 

The plan of the paper is as follows. In section 2 we review some definitions and basic equations 
relegating a description of the numerical scheme we use to solve initial value problems to appendix 
A. Note that no anelastic approximation is made here. In section 3 we go on to describe the results 
we obtained for the response of a simple model of a giant planet as a polytrope of index n = 1 
to tidal encounters. We present results for the total energy and angular momentum transferred to a 
planet without a solid core for the full range of rotation rates and a variety of encounter pericentre 
distances. We then go on to consider planet models with solid cores of radius 25% and 50% of the 
total radius illustrating evidence of the focusing of inertial waves, the emission of shear layers 
at critical latitudes and wave attractors, in these cases. The role of these phenomena has recently 
been emphasised by Ogilvie & Lin (2004) and Ogilvie (2009) for the situation where the planet 
is tidally forced while in a fixed circular orbit. However, in our case with the smaller core, the 
energy and angular momentum transferred hardy differs from what is found for the coreless case. 
In the larger core case the interaction is somewhat stronger with the energy transfer being about 
one order of magnitude larger for the model with a core. 

In section 4 we obtain the normal mode frequencies for the coreless planet model by taking 
Fourier transforms of the time series provided by the post tidal interaction responses and locating 
peaks in the corresponding power spectra. The normal mode frequencies are compared with those 
obtained using the spectral method described in PI and also the WKBJ method described in IPN, 
both of which adopted an anelastic approximation. The spectra obtained from the simulations 
are seen to contain p and / modes in addition to the inertial modes. The frequency locations of 
the main global inertial modes are found to be in good agreement with those obtained by other 
methods for azimuthal mode numbers m = and m = 2. 

In section 4.1.1 we investigate the dependence of the simulation normal mode frequencies 
(or eigenfrequencies) on the planet rotation rate paying a special attention to the case of inertial 
waves. In the anelastic approximation it can be shown that eigenfrequencies of inertial modes are 
proportional to the planet's rotation rate, Vl, so it is natural to express their values in terms of it. This 
proportionality breaks down when the anelastic approximation is not used and the positions of the 
mode eigenfrequencies measured in units of Vl are shifted when it changes. We compare positions 
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of the simulation normal mode frequencies with those obtained using the anelastic spectral method 
but with a correction obtained by applying the first order perturbation theory given in IPN. The 
agreement is very good apart for one case that is apparently affected by an avoided crossing. 

For completeness we consider the tidal response of a coreless model of the type considered by 
Goodman & Lackner(2009) that is in hydrostatic equilibrium under a fixed quadratic gravitational 
potential in section 5. These authors pointed out that a simple analytic solution is available and 
shows that no inertial mode excitation is expected in this case. Our numerical results indicate that 
any residual inertial mode excitation occurring as a result of numerical effects is very small in this 
case as should be expected. 

In section 6 we compare the total energy and angular momentum transferred to a coreless 
polytropic planet with the anelastic spectral results already given in IP. We conclude that these are 
in good agreement and therefore the anelastic approximation is a valid approximation. Thus we do 
not confirm the suggestion of Goodman & Lackner (2009) that these exchanges could be seriously 
overestimated as a consequence of its use. Finally in section 7 we discuss our results. 

2 BASIC DEFINITIONS AND EQUATIONS 

Here we adopt a spherical polar coordinate system (r, 9, 0) and the associated spherical coordinate 
system with origin at the centre of mass of the planet. 

We consider the response to a tidal perturbation with associated gravitational potential "^ext 
written as a Fourier series in terms of the azimuthal angle in the form 

oo 

"^ext = TZeY^ ^ezt,mexp(im0). (1) 

m=0 

Here TZe denotes that the real part is to be taken and although the sum is in general over all m, we 
consider only m = and 2 here. Similar Fourier decompositions are used for the state variables 
of the planet. In the linear approximation each Fourier component responds individually to the 
corresponding component of the tidal potential. 

The planet is characterised by its mass M*, radius R^ and the associated frequency 



where G is the gravity constant. The associated energy and angular momentum scales are E^ = 
GM^/R^ and L* = M^y^GM^R^ respectively. For definiteness, we consider below a planet hav- 
ing approximately Jovian values of mass and radius: M^,=2 x 10^^ g and i?^, = 7 x lO^cm. Accord- 
ingly, it is assumed that E^ ^ 3.81 x iC^, L^ ^ Q x iC^ and i7* ^ 6.23 x 10"^ in cgs units. 
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Note, however, that our resuks, given in dimensional form, may be scaled to other values of M* 
and i?* when expressed in the natural units introduced above. 

Following our previous work (Papaloizou & Ivanov 2005, hereafter PI, Ivanov & Papaloizou 
2007, hereafter IP) we assume in this paper that the planet moves on a highly eccentric (formally, 
parabolic) orbit around a source of gravity of mass M. In the parabolic limit the orbit may be 
characterised by its pericentre distance, Rp. The simulations carried out in this paper started with 
the perturbing mass at a distance eight times the pericentre distance. Tidal forces being propor- 
tional to the inverse cube of the distance are negligible beyond this point (see eg. Faber et al 2005). 
For fixed masses, as an alternative to the pericentre distance, in this paper we define an encounter 
through the parameter 



r^=^M,Rl/{MRl), (3) 

first introduced by Press & Teukolsky 1977, hereafter PT. Note that this expression for r] is strictly 
speaking valid only when M^ <^ M. A more general expression is given in equation (lAlOl ) of 
appendix A. 

2.1 Perturbed equations of motion 

We assume that the planet is rotating with uniform angular velocity fi directed perpendicular to 
the orbital plane. The hydrodynamic equations for the perturbed quantities take the simplest form 
in the rotating frame which we use later on. 

Since the planet is fully convective, the entropy per unit of mass of the planetary gas remains 
approximately the same over the volume of the planet, and the pressure P can be considered as a 
function of density p only, P = P{p). The perturbations of the planet may be considered as adia- 
batic, and therefore the same functional dependence P = P{p) holds during perturbation as well. 
This condition is often referred to as the barotropic condition. In the barotropic approximation the 
linearised Euler equations take the form (see PI) 

-^ + 2rj X V = -SJW + -, (4) 

ot p 

where 

W = clp/p + -^,,t, (5) 

V = (t>r, f e, V(f,) is the Eulerian velocity perturbation, p is the Eulerian density perturbation, c^ = 



dP/dp is the adiabatic sound speed, iy is the viscous or diffusive force pr unit volume and as 
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indicated above \&e^.t is the external forcing tidal potential. The velocity perturbation v is related 

to ^ = (^r, ^e, ^0), the Lagrangian displacement vector, through 

The linearised continuity equation gives 

P' = pcl = -clV-{pi), (7) 

where P' is the Eulerian pressure perturbation. Note that the centrifugal term is absent in equation 
(j4]) being formally incorporated into the potential governing the static equilibrium of the unper- 
turbed star and there is no unperturbed motion in the rotating frame. 

As in our previous work (eg. Papaloizou & Pringle 1981, PI, IP) we shall neglect centrifugal 
distortion of the basic equilibrium which will enable us to adopt a spherically symmetric unper- 
turbed model and density distribution. We additionally simplify matters by adopting the Cowling 
approximation which neglects perturbations to the planet gravitational potential, being equiva- 
lent to regarding the planet as moving in a fixed specified spherical potential. Although we shall 
adopt a standard polytrope with ra = 1 as the basic unperturbed model which is not very centrally 
condensed, the main focus is on the low frequency inertial modes for which self-gravity is not 
expected to play a very important role. This is especially the case for small scale modes and waves 
of the type that arise when a solid core is considered. The Cowling approximation should also be 
valid when estimating the effects of compressibility on these kinds of disturbances. 

Provided that the expressions for the density, sound speed are specified for some unperturbed 
model of the planet as well as the form of the diffusive forces (see below) the set of equations 
(jH — E]) is complete. Details of its numerical solution are given in appendix A. 

2.2 Addition of diffusive effects 

Although problems of the type considered here are linear, because the spectrum is in a general 
sense singular, arbitrarily small scales can develop over time. It is important to note that such 
behaviour has a physical basis and does not arise from pure numerical artifacts. 

For example the boundaries of the simulations we perform consist of the polar axis, the equa- 
torial plane and the surface of the planet. Because of the coordinate singularity at ^ = 0, the polar 
axis is replaced by the cone 6 = On^^ , see below, the latter angle being small. The angle between 
this cone and the equatorial plane is accordingly slightly less than a right angle. As shown by 
Ralston (1973), inertial waves with an appropriate frequency can undergo a sequence of alternat- 
ing reflections at these boundaries, resulting in wave action accumulating in the corner, where the 
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Energy transfer as a function of Omega 




1.5 
time { days) 

Figure 1. The energy measured in tlie rotating frame transferred to the planet as a function of time due to the m = 2 tide is plotted for 
different angular velocities for an encounter with rj = 4. The curves from the uppermost to lowermost at large times are for the non rotating case, 
il/ilt, = 0.16, Cl/ilt = 0.36 ( the curve obtained with twice the numerical resolution 400*400 is superposed and can barely be distinguished), 
fl/flt = 0.51, Cl/flt, = 0.8 and f!/f!» = 1.12 respectively (the latter case was considered for the purposes of illustrating trends). In this and 
other similar Figures below the time unit is days and the energy is expressed in units of the natural energy scale of the planet Et, . 



boundaries intersect. Hence a reflectosingularity is produced implying that such frequencies are in 
the continuous rather than point spectrum. Note that even small perturbations to boundaries can 
introduce effects of this type. 

In order to avoid numerical problems arising from such phenomena we incorporate a diffusivity 
in the simulations which is adequate to effectively eliminate such artifacts in models with no 
specified central cores. 

We added simple diffusive forces per unit volume in spherical coordinates of the form 



d f ^^r\ I d ( . Jy\ 
or \ or J r^ ofi \ Oji I 



(8) 



where the diffusion coefficient or effective kinematic viscosity v = uoy/GM^R^ was taken to be a 
constant. For the simulations reported here with coreless models and a grid resolution of 200 x 200 
(for more information about the numerical grid see appendix A), we adopted z/q = 9.36 x 10^^, 
and for these models with a grid resolution of 400 x 400 we adopted z/q = 4.68 x 10^^. For models 
with core radius 0.25i?* we adopted uq = 4.57 x 10^^ at a grid resolution of 200 x 200 and 
z/q = 2.28 X 10^^ at a grid resolution of 400 x 400. Finally, for models with core radius 0.5 we 
adopted z/q = 1-7 x 10"^ at a grid resolution of 400 x 400 and z/q = 8.5 x 10^® at a grid resolution 
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2.0j»1O~'' 



Angular momentum transfer as a function of Omega 




1.5 
time { days) 



Figure 2. The angular momentum transfen'ed to the planet induced by the m = 2 tides is plotted as a function of time for different angular 
velocities for an encounter with r; = 4. The curves from the lowermost to uppermost at large time are for Sl/Q* = 1.12, f!/C = 0.8, 
n/f2, = 0.51, n/f2, = 0.36 (curve for 400 X 400 resolution superposed), Q/Q, = 0.16 and f2 = respectively. In this and other similar 
Figures below the time unit is days and the angular momentum is expressed in units of the angular momentum scale of the planet L*. 

of 800 X 800. With these choices the diffusivity decreases as the resolution increases. We remark 
that energy and angular momentum transfers measured by evaluating the canonical energy and 
angular momentum (see appendix A) just after the encounter are robust to changes in numerical 
resolution and diffusivity as they should be (see below). 



3 A GENERAL DESCRIPTION OF NUMERICAL RESULTS 

3.1 The energy and angular momentum transfer for a close encounter 

We here give results for the energy and angular momentum transfer as a result of a tidal encounter 
at fixed pericentre distance for different angular velocities. We have considered a range of i^/fi*, 
where Vt^^ is defined in equation ^, up to 1.12 even though we have neglected centrifugal terms 
in the model. This is done in order to more fully explore and gain understanding of the functional 
dependence of the response on its defining parameters. The introduction of centrifugal terms would 
not change the essential response problem although deviations of the form of the effective surface 
boundary may play some role (see eg. discussion in section |Z2] above) . 

A first set of runs was carried out for r^ = 4 taking into account the m = 2 component of 
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Figure 3. Contour plots for the real parts of the m = 2 Fourier components of W (upper left panel), ^r (upper right panel), ^g (lower right panel) 
and §0 (lower right panel) after t = 2.523days for the model with rj = 4 and O/H* = 0.36. In this and other similar Figures, lengths are 
expressed in cgs units. For reference ij, = 7x10^. 



the tidal forcing. The contribution from this has found to dominate that of the m = component 
(see below). In all of these and other similar cases with no interior solid core, after the encounter 
the canonical energy and angular momentum attain either expected constant values or values that 
slowly decay on account of the imposed numerical diffusion. We have considered resolutions of 
200 X 200, 400 X 400 and for models with solid cores 800 x 800. As illustrated below the energy 
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Energy transfer as a function of Omega 




time { days) 



Figure 4. The energy transferred to the planet as a function of time due m = 2 tide is plotted for different angular velocities for an encounter with 
rj = 8\/2- The curves from the uppermost to lowermost at large time are for fl/flt, = 0.36, fl/fl* = 0.27 (the curve obtained for f!/f!» = 0.45 
overlaps), U/Q, = 0.6, n/U, = 0.17, U/U, = 0.8, U/U* = 0.1 and U/U, = 0, respectively. 



Anguiar momentunn transfer 03 o function of Onnega 
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Figure 5. The angular momentum transferred to the planet is plotted as a function of time for different angular velocities for an encounter with 
rj = 8\/2. The curves from the uppermost to lowermost at large time are for f2/f2» = 0.1, f2 = 0, f2/f2, = 0.8, f2/f2, = 0.6, f2/f2, = 0.17, 
n/n» = 0.45 and f2/f2, = 0.27 ( the curve obtained for f2/f2, = 0.36 overlaps), respectively. 
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Figure 6. The energy transferred to the planet as a function of time due to the axisymmetric tide (m = 0) is plotted for different angular velocities 
for an encounter with rj = A. The lowermost curve at large times corresponds to r2/f2, = 0.8, the curve immediately above this to f2/f!» = 0.36 
and curve immediately above that , which is almost coincident, is for the non rotating case. For purposes of comparison the uppermost curve is for 
the non rotating case with the m = 2 component of the tidal potential. 



Energy transfer for core radius 0.25R 
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Figure 7. The energy transferred to the planet as a function of time due to the m = 2 tide is plotted for a model with a solid core of radius 0.25R* . 
The encounter is for r] = 4 and f2/f2, = 0.8. The curves from the lowermost to uppermost at large time are for a resolution of 200 X 200, a 
resolution of 400 * 400, for comparison a model with no core and a model with no core and no rotation. 
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Figure 8. The angular momentum transfen'ed to the planet as a function of time due to the m = 2 tide is plotted for a model with a solid core of 
radius 0.25/?, . The encounter is for r) = A and Q/C = 0.8. The curves from the lowermost to uppermost at large time are for a resolution of 
800 X 800, a resolution of 400 * 400, for comparison a model with no core and a model with no core and no rotation. 



and angular momentum transferred between central object and planet is measured just after the 
encounter. This measures the physical exchange and is numerically converged. 

The energy transferred to the planet as a function of time due to the m = 2 tide is plotted 
for different angular velocities for an encounter with r/ = 4 in Figure \T\ Results are given for 
the non rotating case, fi/fi* = 0.16, 0.36, 0.51, 0.8 and fi/i7* = 1.12, respectively. The angular 
momentum transferred to the planet as a function of time for the cases illustrated in Figure [T] are 
plotted in Figure[2l The energy transferred is seen to decrease monotonically with angular velocity 
by about a factor of thirty as fi varies between zero and the maximum value considered. The energy 
transferred is expected to decrease as the star is spun up from zero (Ivanov & Papaloizou 2004) 
but may increase again at large values of r^ (see PI, IP and below). 

The angular momentum transferred to the planet has a maximum absolute value for the non 
rotating model and decreases with increasing angular velocity, becoming zero when f^/fi* = 0.51. 
For larger values of n angular momentum is transferred from the planet to the star. Such a sign 
change is indeed expected when fi is comparable to the angular velocity at pericentre, see PI and 
IP. We remark that the simulations for fi/f]* = 0.36 were done at resolutions of 200 x 200 and 
400 X 400 and are illustrated in Figures [U and [2l The results of these are indistinguishable. 
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The form of the planet response in these cases is found to be global. To illustrate this, contour 
plots for the real parts of the m = 2 Fourier components of W^ ^r, ^e and ^<^ (lower right panel) 
after t = 2.523days, well after the encounter is complete, for the simulation with fi/f]* = 0.36 
are shown in Figure [3] These have global form that is typically found at any time. As expected, 
there is no evidence of evolution to small scale structures indicative of the emission of shear 
layers at critical latitudes and wave attractors (eg. Ogilvie & Lin 2006, Rieutord & Valdettaro 
2010) in these cases. The global form is indicated both by the accurate conservation of canonical 
energy and the accurate conservation angular momentum even though some numerical diffusion 
was present and also that the energy and angular momentum transferred is in good agreement with 
values obtained using a basis function approach developed in our previous papers PI and IP with a 
relatively small number 15 x 15 of trial functions, see the next section. 

3.2 The energy and angular momentum exchanged for a more distant encounter 

We now consider simulations of the above form but with a larger value of r^ = 8v^. In addition 
to the non rotating case, the angular velocities considered were i^/f]* = 0.36,0.27, 0.45,0.6, 
0. 17, 0.8 and 0.1. The energy transferred to the planet as a function of time is illustrated in Figure 
in and the angular momentum transferred to the planet is illustrated in Figure [5] . In this case the 
energy transferred initially increases with angular velocity but attains a maximum value before 
decreasing again. The reason for this behaviour is that the / mode is ineffective at this value 
of r] leading to a weak interaction for the non rotating model. However, once the star is spun 
up, the interaction can strengthen through the inertial mode response. In this case the maximum 
transfer occurs for fi/f]* ~ 0.36. Similarly the angular momentum transfer is very small for the 
non rotating model It increases to a maximum for Vl/Q^: = 0.1 and then passes through zero at 
r2/i7* ~ 0.13. For larger angular velocities angular momentum is transferred from the planet to 
the star, this transfer being largest in magnitude when fi/i7* ~ 0.32 beyond which it decreases in 
magnitude. The sign change again occurs at a value of the angular velocity comparable to that at 
pericentre passage and the planet response remains global. 

3.3 The axisymmetric case m = for rj = A 

In order to evaluate the effectiveness of the m = tide we carried out simulations using the m = 
component of the tidal potential for r/ = 4. The energy transferred to the planet as a function of 
time is plotted for fi/fi* = 0.8, 0.36 and the non rotating case in Figured The non rotating case for 
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Figure 9. As in Figure[3]but for a model with core radius of 0.25-R, and Sl/Sl, = 0.8 after t = OAOSdays. 



m = 2 is also shown. The energy transfer for m = is always found to be an order of magnitude 
or more smaller than the corresponding m = 2 contribution. Therefore this can be neglected. 
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Figure 10. As in Figure|9]but at t = 0.776days and with contour line plots in order to illustrate the inertial wave propagation pattern 



3.4 Planet models with a solid core 



3.4. 1 A planet with core radius equal to 0.25/?=^ 

We now go on to consider planet models with a solid core. In these cases the long term behaviour 
differs from the coreless models owing to the development of wave patterns related to critical 
latitude phenomena and wave attractors (eg. Ogilvie & Lin 2006, Ogilvie 2009 and references 
therein) that correspond to the development of increasing concentrations of wave energy along 
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Figure 11. The energy transferred to the planet as a function of time due to the m = 2 tide is plotted for a model with a solid core of radius O.S-R*. 
The encounter is for rj = 4 and f2/f2, = 0.8. The curves from the lowermost to uppermost at large time are for a comparison model with no core, 
for a resolution of 400 X 400, for a resolution of 800 * 800, for further comparison a model with no core and no rotation. 



closed characteristic paths. However, it is important to note that our calculations are not for fixed 
frequency forcing, which is the situation considered for the discussion of critical latitude phenom- 
ena and wave attractors. Here we are concerned with a combination of different frequencies that 
should lead to a broader pattern. However, we still might expect that singular behaviour develops 
asymptotically for large time, in the non diffusive case. In the diffusive case this would be limited 
but energy dissipation would occur at an enhanced rate as compared to the non singular coreless 
models. This feature is observed in our simulations. 

The energy transferred to the planet as a function of time due to the m = 2 tide is shown in 
Figure|7]for a model with a solid core of radius 0.25-R* for r] = A and Vt/Vt^: = 0.8. Simulations 
have been performed with resolutions of 200 x 200 and 400 * 400. The results agree closely with 
some small additional dissipation in the lower resolution case, which is to be expected. Interest- 
ingly, the energy transferred to the planet is almost the same in the cases with and without a core. 
This transfer is less in magnitude than in the non rotating case, which gives almost identical results 
to those found for the non rotating coreless case. 

The corresponding angular momenta transferred to the planet are plotted as a function of time 
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Figure 12. The angular momentum transferred to the planet as a function of time due to the m = 2 tide is plotted for a model with a solid core 
of radius 0.5i?, . The encounter is for rj = A and C/C» =0.8. The curves from the lowermost to uppermost at large time are for a resolution 
of 800 X 800, for a resolution of 400 * 400, for the con'esponding model with no core and for further comparison a model with no core and no 
rotation. 



in Figure [8l As for the energy transferred, the angular momentum transferred in this model and the 
coreless model are very similar. 

The global response pattern after t = OAOSdays is illustrated in Figure[9l This is shortly before 
the end of the encounter. Although there is evidence for some small scale structures developing 
near the core boundary, the basic response appears global with little evidence of critical latitude 
phenomena or wave attractors. However, this appears at later times. Response pattern contours are 
shown in Figure [10] at t = 0.77days, being well after the end of the encounter. Localised inertial 
wave propagation and reflection patterns are clearly visible. As expected these present a constant 
angle of propagation to the vertical and graze the core boundary. 

The fact that the development of these structures occurs at late times is indicative of why the 
coreless and cored results agree in this case. We note that this may be expected when the time 
scale for the encounter is not long compared to the rotation period. In such a case there would not 
be adequate time for inertial waves to propagate, reflect and develop structures indicating shear 
layers or attractors. This might be expected to be especially the case when the core is small and so 
presents a small target for propagating waves. Physically, as far as the inertial mode excitation is 
concerned, the encounter would appear impulsive in which case the energy and angular momen- 
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Figure 13. As in Figure [Tol but for tlie model with a core radius of 0.5/?* at t = 0.5S9days. 

turn transfers should not be affected by the later development of shear layers or wave attractors. 
Therefore, as long as the core is not too large we expect the same result as for a coreless model. 
These ideas are developed from a mathematical point of view in appendix B. Although the en- 
counter time scale and rotation period are comparable in this case, in accordance with the above 
discussion, we do not see the development of wave propagation patterns until the late stages of 
the encounter, and the cored and coreless model energy and angular momentum transfers are very 
similar. 
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Figure 14. As in Figure|9]but for a model witli core radius 0.25-R* , f2/f2» = 0.36 and rj = 8\/2 at t = 2.326days, wliicli is after the encounter. 
Contour line plots illustrate the inertial wave propagation pattern 



3.4.2 A planet with core radius equal to 0.5i?* 

We have also performed simulations a model with a solid core of radius 0.5-R*. Again the encounter 
was for 7] = A and Vt/Vt^ = 0.8. The energy and angular momentum transferred as a function of 
time are shown in Figures [H] and [l2l respectively. Results are given for resolutions of 400 x 400 
and 800 x 800. These agree closely but with again more dissipation in the lower resolution case. 
In contrast to the case of the smaller core the energy transferred is larger than in the coreless 
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Figure 15. The real part of the 9— component of velocity in an inner point of the planet as a function of time, rj = 4\/2 and m = 2 for these tidal 
encounters. Simulations for Q/f2, = 0, 0.24 and 0.36 are plotted with dotted, solid and dashed curves, respectively. 



case by a factor of ~ 8 — 9. This is probably because the larger core enabled a wave structure to 
develop faster. To illustrate the wave pattern, contour plots of the response are shown in Figure [T3] 
at t = 0.59days, which is shortly after the encounter. These appear to be of a similar form to, but 
somewhat more intricate than the model with a smaller core. 



3.4.3 Cored models for slower rotation and larger r] 

In order to indicate that the results described above are not radically altered when Vl is reduced 
and T] is increased, we also simulated cored models with r] = 8v^ and fi/fi* = 0.36. The models 
with core radii 0.25-R* and 0.5-R* behaved very similarly, in relation to their coreless counterparts, 
to those described above for larger fi and smaller r]. Figure [M] illustrates the wave propagation 
pattern that results from the encounter when the core radius was 0.25/?=,,, which is indeed similar 
to that found for the cases with larger rotation. 
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Figure 16. Fourier spectra obtained by the Fourier transform of tlie velocity data as functions of tlie ratio a/Clt for rj = 4\/2. Tiie case m = 2 
is sliown. Tlie solid (black) curve corresponds to the non-rotating case £7 = while the dashed (red) and dash-dotted (blue) curves are for 
fl/Qt = 0.24 and 0.6 respectively. The vertical lines delimit the range of a, where inertial modes are theoretically possible, with line styles 
identical to those used for corresponding values of Q. 

4 A COMPARISON OF THE NUMERICAL, SPECTRAL AND WKBJ APPROACHES 

4.1 Normal mode frequencies obtained from the simulations and their comparison with 
the basis function and WKBJ approaches 

Here we investigate the spectrum of normal modes excited in the above simulations determined 
from fast Fourier transform of the time series. We compare the derived inertial mode frequencies 
with the results obtained in IPN. 

As time series used for our Fourier analysis, for definiteness, we consider the time dependence 
of the real and imaginary parts of the ^—component of velocity at an inner point of the planet, with 
r = 4.8 X lO^cm and 9 = 67 degrees, after a tidal encounter with r] = 4v^, for different values 
of the angular velocity Vl. This is shown in Fig. [15] for fi/fi* = 0.24 and 0.36. The sharp change 
occurring around 0.5 day is due to the tidal encounter itself, the pulsation modes are generated 
during this period of time. After approximately 1 day the influence of the perturbing tidal field 
becomes unimportant and the planet's perturbation becomes quasi-periodic. As follows from Fig. 
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Figure 17. The Fourier spectra as functions of ttie ratio cr/Q, with only the range corresponding to the inertial branch of pulsations |o"/n| C 2 
explicitly shown, for m = 2 and r? = 4\/2. The solid (blaclc), dashed (red) , dotted (green) and dot-dashed (blue) curves are for f2/f!» = 
0.24, 0.36, 0.6 and 1.0. Note that the latter case is unrealistic but included for completeness. The symbols show positions of the WKBJ modes 
obtained in IPN, with different style corresponding to different WKBJ mode order, I. The circles, squares, diamonds and triangles are for i = 1, 2, 3 
and 4, respectively. 



[T5l in the case of the non-rotating planet, only one relatively short period of pulsation can be 
clearly seen. This is the period of the fundamental prograde mode. When the planet rotates several 
periods can be seen with one short period corresponding to the fundamental mode together with 
several relatively long periods determined by presence of inertial modes in the spectrum. 

To obtain the Fourier transforms we removed the first part of the data series to eliminate the 
contribution of direct forcing by the tidal field, typically of duration 5 x 10"^ — lO^s. The results 
are presented in Figs. [T6lfT8l In Fig. [16] we compare the spectrum obtained for the non-rotating 
planet with those obtained for the rotating cases, for m = 2. The spectra are plotted as functions 
of the ratio a/fi*, where a is the frequency. The peaks in Figure [T6] correspond to positions of 
normal modes. The black curve represents the results obtained for fi = 0. The positions of the 
two peaks are symmetric about a = 0, though the peak amplitudes are not. This is because the 
prograde modes with positive values of a are excited much more effectively than the retrograde 
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modes, for a non-rotating planet [J. The prominent peaks with a ~ ±1.6^2=,, give the positions of 
the fundamental mode. Note that their eigenfrequencies are larger than those calculated for the 
same planet model but without the Cowling approximation. The peaks with larger values of a /Vt^, 
give the positions of p-modes (pressure waves). Note the absence of any peaks in between the two 
peaks corresponding to the fundamental mode. This is because there is no branch of low frequency 
oscillations in a non-rotating barotropic planet. When the planet rotates the peaks corresponding 
to fundamental modes are shifted leftwards, with the magnitude of the shift proportional to Vt. This 
effect is briefly discussed below. 

The most notable difference between the non-rotating case and the rotating ones is the presence 
of peaks in between the two peaks corresponding to the fundamental modes in the rotating cases. 
These are the peaks corresponding to the inertial modes. Note that all these peaks lie within the 
theoretical range allowed for inertial waves \a ^ 2Vt\. This range is marked by vertical lines with 
styles corresponding to those associated with the corresponding angular velocity. For small values 
of VL, significantly less than the largest value Vl/VL^ = 0.6 shown in Fig. [16], there are two notably 
large peaks within the inertial range, one with a positive and one with a negative value of a. These 
are identified with two main global modes discussed in detail in PI, IP and IPN. From the point 
of view of the WKBJ theory developed in IPN, these modes correspond to smallest value of an 
integer / classifying WKBJ modes, called the WKBJ order later on. In the case of Vl/Vl^ = 0.6, 
represented by the dot-dashed curve, there is a third noticeable peak close to a ~ 2fi and within 
the inertial mode range. It is the prograde fundamental mode, which has been shifted within this 
frequency range by rotation. A number of peaks with smaller amplitudes give the positions of 
higher order inertial modes. 

The eigenfrequencies of inertial modes are naturally expressed in units of il and we show the 
part of the spectrum within the inertial range as a function of the ratio a/fi, for the case m = 2 
in Fig. [l7] and for the case m = in Fig. [TSl Additionally, we show positions of WKBJ modes 
calculated in IPN by symbols with larger values of symbol ordinates and sizes corresponding to 
smaller values of the WKBJ order /. As was discussed in IPN among the WKBJ modes with the 
same value of / there are modes, which were not identified with results obtained by other methods. 
These modes are situated close to boundaries of the allowed region a = ±2fi as well as close 
to the origin a = 0. The theory developed in IPN is likely to be beyond limit of its applicability 

^ Let us recall that it is assumed in this paper that that the direction of planet's rotation and the orbital motion coincide and that planet's perturbations 
are proportional to e^H""'"""^), where </< is the azimuthal angle in the rotating frame. It is easy to see that the modes with cr > (ct < 0) are 
prograde (retrograde) with respected to the orbital motion. 
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Figure 18. Same as Fig. ll7l but for the case m = 0, rj = 4. Since in this case the spectra are symmetric about the origin a = only the positive 
values of a within the inertial range cr/Q ^ 2 are shown. Q/Qt = 0.36 for the solid (black) curve and Q/C = 0.8 for the dashed (red) one. 
Symbols show the positions of WKBJ modes. Note that the WKBJ mode order I runs starting from Imin = in this case. The black circles, red 
squares, green diamonds, blue triangles and magenta stars are for i = 0, 1, 2, 3 and 4, respectively. 



close to the boundaries and to the origin a = 0. We expect that these unidentified modes could be 
absent in a more advanced WKBJ theory. Therefore, we do not show positions of modes having 
most negative, most positive and the smallest values of a, for a given value of /. 

In Fig.[l7]the case of m = 2 is shown. The solid black and dashed red curves correspond to the 
cases of relatively small fi/i7* = 0.24 and 0.36, which should be compared with WKBJ results 
obtained in the anelastic approximation. Two global modes are clearly seen with their positions in 
a close agreement with positions of the WKBJ modes having the smallest value of Imin = 1 for 
m = 2. The WKBJ modes corresponding to / = 2 and 3 are also shown. Their positions are close 
to some numerically obtained peaks, notably in the case of the / = 2 modes having the ratio a/fi 
close to —1.5, —0.8 and 1.25. The fact that some of the next order WKBJ modes are not seen in 
the numerical data is not worrying since the amplitude of modes excited after the tidal encounter 
depends of the shape of the Fourier transform of the dependence of the tidal potential on time, 
which, in its turn, depends very sensitively on the value of Q, see PI and IP. With increase of n the 
retrograde modes with negative values of a are excited more efficiently and their amplitude grows. 
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Figure 19. Eigenfrequencies a/Q of a few 'global' inertial modes as functions of the ratio Q/Qp, where flp is is atypical frequency of periastron 
passage, see equation )A9t . The dot-dashed curve having a/Cl ai 1.3 when H — > is the analytical curve for the axisymmetric m = global 
1 = 2 mode. The curves having a/Q Ri 0.5 and —1.1 when f2 — >■ are for the global I = 1 m = 2 modes, while the curve with cr/U ai 0.4 is 
for the next order I = 2 curve. Symbols joined by dotted lines represent the corresponding results of numerical calculations. One can see that the 
agreement between the numerical and analytical results for the global m = mode and the global rn = 2 mode with negative values of cr is quite 
good. The analytical results for the global rn = 2 mode with positive values of a agree with the numerical ones only when fl is sufficiently small. 
In the opposite case of large values of Q the numerical results are close to the analytical curve for the next order m = 2 I = 2 mode. This result 
may be explained by the effect of avoided crossing. 



see Figure [TTl In the case of the largest considered value of f2 = fi^, the peaks corresponding to the 
retrograde mode with / = 1 and 2 are strongly amplified, and the peaks close to cr/fi = — 1.7 and 
—0.7 may be identified with the WKBJ modes having 1 = 3. Also, the peaks are shifted towards 
the origin a = as fi increases, an effect which may be explained by corrections to the anelastic 
approximation, see IPN and below. In the cases of fi/fi* = 0.6 (the dotted green curve) and 
ri/n^: = 1 (the dot-dashed blue curve) there is an additional noticeable peak with a positive value 
of a. This corresponds to the prograde fundamental mode being shifted into the inertial range due 
to rotation. 

In Fig.[T8]we show the case m = 0. Since in this case the spectra are symmetric about a = 
only positive frequencies within the inertial range cr < 2fi are considered. Similar to the previous 
case symbols show positions of the WKBJ eigenmodes calculated by IPN. As in Fig. [16] larger 



ordinates correspond to smaller values of the WKBJ order / with the smallest / = /„ 
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Figure 20. Eigenfrequencies ct/H, of m = 2 fundamental modes as functions of the ratio Sl/f2». The solid and dashed curves represent 
the theoretical results for the prograde and retrograde mode, respectively. Symbols joined by the dotted lines represent the results of numerical 
calculation. One can see that the agreement between the theoretical and numerical results is much better for the prograde mode. Note that the same 
phenomenon was observed in numerical calculations of normal modes of rotating polytropic stars, which do not use the Cowling approximation, ( 
see eg. Managan 1986). 



denoted by circles and the eigenfrequencies of modes with largest shown / = 4 denoted by stars. 
We see from Fig[T8]that the peak with largest amplitude has a value of a close to the identified 
main global mode with / = and a ~ 1.3VI. Two identified WKBJ modes at the next order, also 
have well pronounced peaks with positions close to the values of their eigenfrequencies. One can 
also possibly identify the eigenfrequencies corresponding higher order modes with some peaks 
with smaller amplitude, especially in the region 0.5 < a/i7 < 1, where there is seemingly a good 
agreement between positions of the peaks and eigenfrequencies of modes having / ^ 6. When 
rotation of the planet increases the positions of the peaks are shifted to smaller values of a. This 
agrees with results obtained in IPN, see also below. 
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4.1.1 Change of eigenfrequencies with Vt 

In the anelastic approximation which is frequently used for calculations of eigenspectra of modes 
belonging to the inertial branch, the mode frequencies are proportional to Vt. This approximation 
is not, however, exact, and it discards certain terms proportional to (i7/i7*)^ in the set of equations 
describing pulsations, see eg. IPN for discussion. In IPN a correction to the mode eigenfrequencies 
due to terms not accounted for in the anelastic approximation was calculated. It was shown that in 
the linear approximation in the small parameter {Vl/VL^Y and neglecting a possible mode mixing, 
the eigenfrequencies have the form 

a, = < + «,(fiM)2fi, (9) 

where a^ is the value of a certain eigenfrequency calculated in the anelastic approximation and an 
explicit expression for dimensionless parameter «« can be found in IPN. We calculate a^ for three 
modes corresponding to the case m = 2 using results obtained in PI and IP. For the main retrograde 
global mode having / = /mm = 1, from the viewpoint of the WKBJ method, the eigenfrequency 
cr" ~ —l.lVl and ai ~ 0.15. For the prograde / = 1 main global mode we have 0-2 ~ 0.54O 
and ^2 ~ —0.5, and for the next order / = 2 mode which has a close value of the eigenfrequency 
(Ta ~ 0.41^2 we get 0:3 ^ —0.022. Additionally, we calculate the correction for the m = 2 main 
global mode with / = Imm = 2, which has ai ^ 1.35^ and ai ~ —0.12. Then, the value of at 
given by equation dH) is compared with what is obtained in the numerical calculations. The results 
of this comparison are presented in Figure \T9\ Note that for the curves corresponding to m = 2 
the parameter i] = 4-\/2 while in the case m = 0, r^ = 4, and we recall that 77 = fi*/fip, where Vlp 
is is a typical frequency of periastron passage explicitly defined in equation (IA9I ) of appendix A. 
As seen from Fig{T9] there is good agreement between the numerical results and those calculated 
with help of equation ^ for the m = global mode and the retrograde m = 2 global mode. 
In the case of the prograde m = 2 global mode only two numerical points with n/Vlp < 1.5 
are in agreement with the analytical result. For larger rotation rates the numerical curve deviates 
from the analytical one calculated for the global mode. Instead, it approaches the analytical curve 
calculated for the next order / = 2 mode in the limit of large rotation. Since the theoretical curves 
corresponding to the global mode and the / = 2 mode intersect each other at n/Vlp ^ 3, the fact 
that at larger rotation rate the numerical results are closer to the curve describing the / = 2 mode 
may be explained by the phenomenon of avoided crossing. The global and the next order mode 
lose their identity at these rotation rates, being in mixed states with mode mixing provided by an 
operator describing the correction to the anelastic approximation, see IPN for its explicit form. 
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Now let us shortly discuss the behaviour of the fundamental modes with change of rotation 
rate. In order to find the corresponding analytical estimate of change of eigenfrequencies we should 
remember that we work in the rotating frame while the fundamental modes are normally studied 
in the inertial frame. In the latter frame there is a simple expression for a correction of the mode 
eigenfrequency, aj, due to rotation 

af - aj = m/Sn, (10) 

where aY is the frequency of a fundamental mode calculated for a non-rotating planet, and (3 is 
determined by some integral over the planet's volume, see eg. Cristensen-Dalsgaard and references 
therein. As was shown by eq. Ivanov & Papaloizou 2004 (3 ^ 0.5 for planetary models, which are 
close to the standard n = 1 polytrope. Eigenfrequencies associated with the rotating frame are 
shifted with respect to the ones in the inertial frame according to the rule cr — t- a — mVl. Therefore, 
in the rotating frame, for the case of the m = 2 mode, we obtain a simple relation 

af = aY-Q, (11) 

where we set /3 = 0.5. This relation is compared with the numerical results in Fig. [201 One can 
see from Figj20] that the numerical and analytical approaches are in a good agreement for the 
prograde mode while in the case of retrograde mode the results differ. The fact that the simple 
approximation (fTTI) doesn't work well for the retrograde fundamental modes, for finite values of 
rotation rates, has been known for quite a long time, see eg. Managan (1986). Since this issue is 
not directly relevant to the purposes of this paper we do not discuss it in more detail here. 

5 CORELESS MODELS UNDER A FIXED QUADRATIC POTENTIAL 

Goodman & Lackner (2009) considered special barotropic models which are in hydrostatic equi- 
librium under a fixed quadratic gravitational potential \E' in the unperturbed state. Thus 

1 dP _ cj. dp _ d^ 

p dr p dr dr ' 

with ^ = ci + C2r^, with Ci and C2 being constant. Such models have a special property in common 

with the uniform density model that {l/zu)d'^ /dzu = {1 / r)d'^ / dr = C2 = uJq is constant. Such 
spherically symmetric models may be found for an arbitrary specified density distribution if the 
pressure is allowed to be determined a posteriori by the condition for hydrostatic equilibrium. 

Goodman & Lackner (2009) showed that these models had the property that when full com- 
pressibility was retained, no inertial modes are excited by a quadrupole m = 2 tidal potential. 
Thus for these models the response is determined entirely by the / modes. In appendix O we 
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0.24 


1.1940 


1.1965 


-1.6750 


-1.6751 


0.36 


1.0988 


1.0990 


-1.8202 


-1.8222 



Table 1. A comparison between / mode angulai' frequencies obtained from numerical simulations with analytically determined values. The first 
column gives the angular velocity, the second and third columns give the analytic and numerical prograde / mode angular frequencies respectively. 
The fourth and fifth columns give the coiresponding quantities for the retrograde / mode. 

review the compressible case and show that no inertial modes are excited when the anelastic ap- 
proximation we use is employed. Thus no spurious inertial mode excitation occurs as a result of its 
use for these models. Nonetheless it is important that the numerical schemes we employ correctly 
represent the tidal excitation of these models. This we demonstrate below for the finite difference 
approach adopted in this paper. The situation with respect to the spectral approach adopted in 
previous papers, which is also found to behave correctly is discussed in appendix O 

5.1 Simulations with quadratic potential models 

The models that we have adopted in order to compare their tidal response to that of the standard 
polytropic model considered in this paper have the same mass and radius as but a density distribu- 
tion that depends quadratically on radius. The central density is given by 

This is smaller than the corresponding value for the polytrope in equilibrium under a self-consistent 
gravitational potential given in appendix |A] because the quadratic potential model is less centrally 
condensed. The central pressure is given by 

IhGMl 

The density and pressure in the interior are then respectively given by p = Pc(l — r^/_R^)and 
P = Pc(l — r'^/Rt)'^- Thus this model also satisfies the equation of state for a polytrope with 
index n = I. However, hydrostatic equilibrium confirms that this is a quadratic potential model 
with potential determined to within an arbitrary constant by \E' = GM^r"^ /{2R^). 

We have performed simulations of quadratic potential models for r] = 4v^ for a number of 
rotation frequencies and a resolution 200 x 200 that can be directly compared to those for the 
corresponding standard polytropes with n = 1 (the global simulation parameters for the different 
models are the same). In particular, we establish that although there is some mild distortion due 
to the form of the computational grid and numerical truncation, that causes some weak inertial 
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Figure 21. This figure illustrates the behaviour of the quadratic potential models discussed in the text with some comparisons to a corresponding 
standard polytropic model for encounters with rj = 4\/2. The left uppermost panel gives contours of the real part of W and the right uppermost 
panel contours of the real part of rvg 13.38 days after pericentre passage for the quadratic potential model with f2/n* = 0.36. The central panels 
give the real part of j^g evaluated at r = 4.8 X lO'-* cm and S = 67 degrees as a function of time between six and eight days after pericentre passage. 
The left central panel illustrates the standard polytropic model (solid curve) and the corresponding quadratic potential model (dashed curve) for 
Q/Q* = 0.24. The right central panel illustrates the same quantities for f!/f!» = 0.36. The lowermost left panel shows the energy transferred 
from the orbit relative to its value for Q = as a function of Q in sec^^ . We recall that f2* = 6.23 X 10^^ sec ~^ . The standard poly trope is 
illustrated with the solid curve and the quadratic potential model with the dotted curve. The lowermost right panel gives the con'esponding plot for 
the angular momentum transferred from the orbit. 
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Figure 22. This shows a comparison of the Fourier spectra obtained by the Fourier transform of the velocity data as functions of the ratio cr/il, for 
r} = 4\/2 for standard polytropic and con'esponding quadratic potential models. The case m = 2 is shown. The solid (black) curve corresponds to 
the quadratic potential model with Q/Q, = 0.24 while the dashed (red) curve is for the coiresponding polytropic model. The dash-dotted (green) 
curve is for the quadratic potential model with f2/f2, = 0.36 and the dotted (blue) curve is for the corresponding polytropic model. The vertical 
lines delimit the range of a, where inertial modes are theoretically possible, with line styles identical to those used for corresponding values of Q. 



mode excitation, the response of the quadratic potential models is determined by the / modes as 
expected. Some numerical results for these simulations are shown in Fig. |2ll 

The lowermost left panel shows the energy transferred from the orbit as a function of il rela- 
tive to its value for i7 = 0. The lowermost right panel gives the corresponding plot for the angular 
momentum transferred from the orbit. We remark that when f] = 0, the energy and angular mo- 
mentum transfers for the quadratic potential models exceed those for the polytropic models by 
factors of 9.2 and 11.0 respectively, because the polytropic models are more centrally condensed. 
The lack of inertial mode response is apparent in the quadratic potential models which do not show 
the related increase in energy transfer for intermediate values of il. The interaction is much weaker 
for intermediate and larger values of f2 for these models. This is also indicated when one considers 
the value of il for which the angular momentum transfer is zero. This defines the angular velocity 
fips for which pseudo synchronisation is achieved (Ivanov & Papaloizou 2004, IP). 

The angular velocity for pseudo synchronisation taking account only / modes was estimated 
in general by Ivanov & Papaloizou (2004) to be given by 
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T] 

However, when inertial modes are taken into account IP find that i^ps ~ l.SSfi^/r/. Thus, assuming 
/ modes dominate the quadratic potential models and inertial modes dominate the polytropic 
models, for the simulations illustrated in Fig. |2T1 we expect fip^ = 071^2* in the former case 
and ilps = 0.27il* in the latter. For the numerical results plotted in Fig. [211 ^ps ~ 0.6fi* and 
fips ~ 0.29^2^, for these quantities respectively. As the estimates of Ivanov & Papaloizou (2004) 
and IP assumed fi to be small, this agreement is satisfactory and implies consistency with the view 
that the polytropic models are inertial mode dominated and the quadratic potential models / mode 
dominated. 

The form of the tidally excited disturbance in the quadratic potential model with fi/f2* = 0.36 
13.38 days after pericentre is illustrated in the uppermost panels of Fig. [211 The same form is 
seen at other times long enough after the tidal encounter. Contours of the real parts of W and rvg 
are shown. The former are approximately vertical as expected from the analytic form W oc w'^ 
(see appendix [O) and the latter are approximately hyperbolae corresponding to the expected form 
rvg oc zuz. 

The central panels of Fig. [2Ilgive the real part of vg as a function of time at an inner point of the 
planet, with r = 4.8 x lO^cm and 9 = 67 degrees, between six and eight days after pericentre pas- 
sage both for the standard polytropic model and the corresponding quadratic potential model for 
fi/fi* = 0.24 and fi/fi* = 0.36. It will be seen that the quadratic potential model responses show 
high frequency behaviour corresponding to excitation of the / modes whereas the corresponding 
polytropic models show in addition a dominant contribution from longer period disturbances that 
are associated with inertial modes. 

In order to better quantify this we have performed a Fourier analysis of longer spans of these 
time series in order to determine the mode spectra for these cases. The results are plotted in Fig. 
[22l From this figure it is apparent that when corresponding models are compared the inertial mode 
amplitudes are an order of magnitude or more smaller for the quadratic potential models, corre- 
sponding to a mode energy a hundred times or more smaller. For the standard polytropic cases 
plotted, the response is dominated by the inertial modes in contrast to the / mode dominated 
quadratic potential models. 

As an additional test we compare the numerically determined / mode frequencies for the 
quadratic potential models to the analytically determined values. For the quadratic potential mod- 



els the / mode angular frequencies are given by a = —Q, ± \/2ujq + ^2^ (Goodman & Lackner 
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Figure 23. The dimensionless quantities e and A determining tiie energy transfer in tlie rotating frame and the angular momentum transfer, 
respectively, are plotted as functions of f2/f2p. The case of a tidal encounter having r; = 4v'2 is shown. The uppermost solid, dashed and 
dotted curves describe results of theoretical calculation of e while the corresponding lowermost curves determine theoretical values of A. Symbols 
represent the results of numerical calculations. The solid curves are calculated with all inertial and fundamental modes accounted for, the frequency 
connections to the fundamental modes given by equation jilt and to two main global modes given by equation J9) are taken into account as well. 
The dashed curve differ from the solid ones by the fact that only the two main global modes are taken into account in the expressions determining 
the contribution of the inertial modes. The dotted curves show the case when the contribution of the fundamental modes as well as the frequency 
corrections are neglected and only two main global inertial modes are taken into account. Note that the dotted curves do not depend on value of r]. 
These curves are, therefore, called 'the universal curves' later on. 



2009). Here the positive sign corresponds to the prograde mode and the negative sign to the ret- 
rograde mode. We compare these values to the numerical values we obtained from the spectral 
analysis of the time series indicated above for i7 = 0, n/n^ = 0.24 and Vl/Vl^ = 0.36 in table 
m The numerical and analytic values correspond to within a relative discrepancy of 2 x 10^^ and 
so are in good agreement. We conclude that our numerical treatment adequately represents the / 
modes and their dominance of the tidal response for the quadratic potential models. 
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Figure 24. Same as Fig. 123 1 but for a closer encounter with r] = 4. For simplicity only two mail global modes are taken into account in the 
expressions determining the energy and angular momentum transfer due to inertial waves. Both contributions of the fundamental modes and of 
inertial modes as well as the corresponding frequency corrections are used to calculate the solid curves while the dashed curves are the universal 
ones. 



6 COMPARISON OF THE ENERGY AND ANGULAR MOMENTUM TRANSFER 
WITH THE RESULTS OBTAINED BY THE SPECTRAL METHODS 



6.1 Numerical simulations versus the spectral approach 

Let us compare the energy and angular momentum transferred for individual simulations with 
results obtained from the basis function (or spectral) approach developed in PI and IP. The ex- 
pressions for energy and angular momentum transfer determined by the inertial modes, AEm and 
AL2 obtained by PI in the anelastic approximation are for convenience reproduced in appendix B, 
see equations (IBll) and (|B2I) . Let us recall that m is the azimuthal mode number. As we discussed 
above the numerical results show that modes with m = give a negligible contribution to the 
energy transfer, and, therefore, we set m = 2 later on. Also note that the quantity AE2 is defined 
in the rotating frame. The energy transfer in the inertial frame is given as AEj = AE2 + f2AL2. 
Contrary to AE2, AEi can be negative. 

It follows from equations (IBll) and (|B2I) that dimensionless quantities 
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Figure 25. Same as Fig. l24l but for a more distant encounter with rj = 8v^- 



e = ri''^E2/E, X = r/^ALa/L, 



(16) 



depend only on the dimensionless ratio Vl/Vlp in the anelastic approximation. Therefore, they con- 
veniently represent the energy and momentum transferred and we compare these quantities with 
what is given by numerical calculations below. In our spectral approach we use the results obtained 
by IP to calculate the overlap integrals (see equation (IB3I) for the definition) and eigenfrequencies, 
for their model of a cold Jupiter with smoothed density jump in a place of phase transition between 
molecular and metallic hydrogen, as this model is like a standard polytrope with index n = I. 
Formally, the quantities e and A are given by infinite summation series over all eigenmodes with 
m = 2. In practice, the number of modes is limited by the resolution of the numerical method used 
, which in turn is determined by the number of basis functions in the case of a spectral method. 
As we mentioned above, in the calculations of PI and IP 225 basic functions were used giving 
450 eigenfrequencies. However, as was mentioned by PI and IP only two main global modes (a 
prograde one with a ^ 0.5i7 and a retrograde one with a ~ — fi give the main contribution to 
the transfer of energy and angular momentum, for small and moderate ratios of Q/Qp. Thus, to 
compare the results obtained by different methods we consider two approaches: the first takes into 
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account all of the eigenmodes available in the spectral method, the second takes into account only 
the contributions of the two main global modes. 

In order to compare the spectral and numerical results we should add to the expressions (fT6l) 
the corresponding contributions arising from excitation of the fundamental modes, for a rotating 
planet. These have been derived in Ivanov & Papaloizou 2004 (see also Lai 199713- Additionally, 
we take into account the frequency correction Q. After these modifications are made, the expres- 
sions (fT6l) acquire dependence on t]. To emphasise the difference between the expressions with 
and without these modifications we call the latter 'universal curves'. 

The results of the comparison of e and A are presented in Figures [23ll25l In Fig. [23] we show the 
case of a tidal encounter with r] = 4-\/2. The solid and dashed curves represent the results given 
by the spectral method with the contribution of the fundamental modes as well as the frequency 
corrections included while the dotted curves are the universal ones. The symbols show the results 
of numerical calculations. One can see that the agreement between the spectral and analytical 
approaches is quite good apart, possibly a region with Cl/Vlp ^ 1.5. In particular, it is clear from 
Figl23]that for a sufficiently small ratio Vl/rip < 1.5 the energy and angular momentum gained by 
the planet are determined by excitation of the fundamental modes while in the opposite case the 
inertial waves dominate the tidal response. Therefore, eg. the curve giving the energy transferred 
to the planet has non-monotonic shape. Initially it decreases with n/Vlp due to suppression of 
the contribution determined by the fundamental modes caused by increase of rotation, then it 
starts to increase due to increasingly more efficient excitation of the inertial waves attaining a 
maximum at il/Vlp ~ 3.5, then it decreases again. One can also see that the role played by the 
frequency corrections to the inertial modes is essential with the universal curves giving much worse 
agreement with the numerical results in the region, where the inertial waves dominate. On the other 
hand it is difficult to see whether the curve with all inertial eigenmodes taken into account (the 
solid curve) or the curve determined by two main global modes (the dashed curve) gives a better 
approximation to the numerical data. We expect that in order to differentiate between these two 
curves even larger rotation rates should be considered. This issue is briefly discussed below, where 
we compare results of the spectral and WKBJ methods. In the rest of this section we consider only 
two main global modes in the expressions for the tidal transfer due to inertial waves. 

In Fig. (|24|) the case of a stronger tidal encounter with smaller value of r/ = 4 is considered. As 



^ Note that the final expressions of Ivanov & Papaloizou 2004 should be multiplied by tt^ due to a misprint in the text, the correct expressions are 
given in an appendix of IP. 
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in the previous case the agreement is very good. As expected the role of the fundamental modes 
is more important at small rotation rate while the role of the frequency corrections becomes quite 
important at large values of rotation. 

Fig. (l25l) shows the case of relatively weak tidal encounter with r] = 8-\/2. In this case the 
inertial waves dominates the tidal response for all non-zero rotation rates considered. Also, the 
role played by the frequency corrections is minor relative to the previous case, and, therefore, the 
curves taking these ones into account are close to the universal curves in this case. Note that the 
numerically obtained value of e is larger than the theoretical value for the maximal Vl/Vlp ^ 9. 
This might possibly be explained by contributions of higher order WKBJ modes to the energy 
transfer, which are not taken into account in the theoretical curve. 



7 DISCUSSION 

In this paper we extended our previous work (IP) on the tidal interaction of a planet on a highly 
eccentric or parabolic orbit about a central star by relaxing the anelastic approximation and con- 
sidering the tidal response of planet models with solid cores. These changes were accomplished by 
carrying out numerical simulations of tidal encounters considered as initial value problems. The 
response of / and p modes was found in addition to that due to low frequency inertial modes. 

We calculated the energy and angular momentum exchanged after pericentre passage for a 
variety of pericentre distances and angular velocities of rotation of the planet. Models both with 
and without a solid core were simulated. In the latter case the state variables showed evidence of 
the emission of shear layers and/or wave attractors after the encounter. 

In the coreless case we studied the spectrum of excited inertial modes and compared it with the 
spectra obtained in IP and with the analytic theory described in IPN. We found good agreement 
for slow rotation. We were also able to account for the variation of the normal mode frequencies 
with the planet's angular velocity. In the case of inertial waves a non-trivial variation is induced by 
terms neglected in the anelastic approximation and it is found to be small for rotation frequencies 
significantly less than critical. Thus our analysis indicated that this approximation is generally a 
valid one for the inertial modes. 

We presented results for the total energy and angular momentum transferred to a planet without 
a solid core for the full range of rotation rates and a variety of encounter pericentre distances 
and found good agreement in the coreless case with results obtained by IP which were obtained 
using an anelastic approximation. When a solid core is introduced the energy exchanged between 
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star and planet is in general not changed by a very large amount. We give both a physical and 
mathematical explanation of this result in appendix B. 

It is a consequence of the good agreement of our results with those obtained previously using 
other methods that the anelastic approximation is valid for the global inertial modes in the slow 
rotation limit. Thus we do not see the large enhancement of the tidal interaction that Goodman & 
Lackner (2009) indicated would occur as a consequence of adoption of the anelastic approxima- 
tion. 

Since in this limit our numerical results fully agree with the previous ones obtained in Ivanov 
& Papaloizou 2004 and IP, where estimates of the circularisation time scales of a planet on a 
highly eccentric orbit have been made, we can state that these estimates are now confirmed by 
direct numerical simulations that do not use the anelastic approximation. 
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APPENDIX A: NUMERICAL METHODS AND SET UP 
Al Computational grid 

The simulations were performed on a n^. x n^ grid (vi, 9j). Here n.,-- and uq are the number of grid- 
points in the radial and theta directions respectively. After numerical experimentation we adopted 
the following set up which was found to give reasonable results for the particular models adopted 
that had an equation of state corresponding to a polytrope with n = I. 

In setting up the radial grid for coreless models we remark that setting the innermost radial 
grid point too close to the origin and/or the outermost point too close to a zero pressure surface at 
r = R^ results in a singular unstable behaviour. To avoid this we set the innermost grid point to be 
ri = 3.77 X 10^^/?^ and the outermost point to be r2 = 0.975-R* and applied the relevant boundary 
conditions there. This choice of innermost point was found to result in regular behaviour near the 
origin while avoiding the real physically singular wave propagation behaviour manifested by the 
the models with sizable cores. Similarly smooth behaviour was maintained close to the bounding 
planetary surface. For models with solid cores ri was replaced by the core radius. The radial grid 
points were then distributed according to 

n = {rf + {i - l)h2Y/^ I = 1,2,. .Mr. (Al) 

Here h2 = (r| — r\)/{nr — 1). 

The angular grid was defined over < ^ < 7r/2 in terms of /Xj = cosOj according to the 
prescription 

/i, = ((j - l)/i^)2/3, j = l,2,...n,,, (A2) 

where h^ = 2/{2nth — !)• The domain 7r/2 < ^ < tt by assuming symmetry with respect to 
reflection of the density distribution in the equatorial plane. 
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Al.l Specification of the state variables on the grid 

The numerical scheme solves for the Fourier components of each state variable defined analo- 
gously to equation ©. These are associated with a particular azimuthal mode number, m, and 
are complex quantities. This is taken as read throughout. The specification of the state variables 
is then staggered so that the density and pressure perturbations p' and P' as well as f^ and ,^(^ are 
evaluated at cell centres {i + 1/2, j + 1/2). The velocity Vr and displacement ^^ are evaluated at 
(i.j + 1/2), and vg and ^g are evaluated at {i + 1/2, j). The specification of the state variables is 
also staggered in time. Thus, ^, P' and p' were evaluated at time levels k,k = 1,2... and v, was 
evaluated at the intermediate levels k + l/2,k = 1,2... 

Al.l Numerical solution 

Equation (HI) is solved by use of operator splitting by writing it in the form 

— = Li + L2 + L3, (A3) 

where Li = —WW , L2 = — 2f2 x v, and L3 = fj,/p and then solving the equations 

^ = L„ z = l,2,3 (A4) 

in sequence. For the cases i = 1 dealing with the pressure forces and external potential a sec- 
ond order time symmetric explicit finite difference scheme was used. For i = 3 dealing with non 
conservative forces, as these in general are of the order of the truncation error, a first order (in 
time) second order in space finite difference scheme was used. For descriptions of similar numer- 
ical schemes used to solve linearised equations see Papaloizou & Pringle (1984, 1987) and Lin 
, Papaloizou & Savonije (1990). The case i = 2 above involving Coriolis forces was integrated 
as a straightforward system of simultaneous equations using a second order Runge-Kutta method. 
Having advanced v, ^ was obtained from equation Q and then to complete the cycle P', and p' 
are found from equation (|7]). 

A 1.3 Boundary conditions 

Boundary conditions need to be applied at the innermost and outermost radial grid points, on the 
equatorial plane 9 = n/2 and on the cone 9 = 9ng . On the inner radial boundary we adopt a 
regularity condition or Vr = 0. On the outer radial boundary we take the Lagrangian pressure 
perturbation to be zero. Finally at both 9 = 9ng and 9 = 7r/2we adopt vg = 0. 
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A1.4 Initial density 

The initial density distribution is taken to be that of a standard polytrope with index n = 1 in 
equilibrium under its self consistent gravitational potential. Thus, for all models, outside of any 
core we adopt 

_ pcR* sin(7rr/-R^) 
P — (Aj; 

Tcr 

where the central density is given by Pc = tcM^/{ARI) and M^ is the total planet mass. 

A2 The form of the tidal potential for m = 2 and m = 

The m = 2 component of tidal potential is 

^e.t,2 = -4(^^§-'^i^'(^) -M-2^<I>,), (A6) 

(see equation ^ where M be the central perturbing mass and Rp be the pericentre distance. The 
corresponding m = axisymmetric component of the tidal potential is 

il„M R„ o , n 

The radial and angular coordinates of the planet are given by 

R = Rp{l + 4 sinh^ x), and cos(0p) = 2Rp/R - 1 (A8) 

respectively. The quantity x is related to the time t through sinh(3x) = 3Q,p{t — to)/lQ, where 



fip = ^G'(M, + M)/i?3 (A9) 

and pericentre passage occurs at t = tg- The simulations carried out in this paper started with the 
perturbing mass at a distance eight times the pericentre distance. Tidal forces being proportional 
to the inverse cube of the distance are negligible beyond this point (see eg. Faber et al 2005). 
For fixed masses, as an alternative to the pericentre distance, in this paper we define an encounter 
through the parameter 



r] = ^M,Rl/{{M + M,)Rl) ~ ^M,Rll{MRl). (AlO) 

The tidally disturbed body is explicitly taken to be a planet of one Jupiter mass with R^ = 7 x 

lO^cm. 

A3 Energy and angular momentum transfer 

The energy and angular momentum transferred during an encounter are found by evaluating the 
canonical energy and angular momentum. These are conserved and well defined in a non dissi- 
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pative system (Friedman & Schutz 1978). In a weakly dissipative system they decay slowly with 
time but their values just after the encounter represent the associated transferred quantities. 

A3. 1 The canonical energy 

When the Lagrangian pressure perturbation vanishes at the outer boundary, the canonical energy 
appropriate to the Fourier mode with azimuthal mode number, m, may be written as 

E, = 0.25(1 + 6m,o) [/ P (ivH + |P'r/(7^)) ^^ - ^ \^r\\dP/dr)ds] , (All) 

where 5mfi) is the Kronecker 5. Here the first volume integral is taken over the planet volume V 
and the second surface integral is taken over the surface area A. The contribution of the surface 
term formally vanishes when the density at the surface is zero. When the density at the surface 
is relatively small as in our numerical model this term is found to give a negligible contribution. 
When tidal forcing operates, the time rate of Ec gives the rate of energy uptake by the planet. After 
the encounter tidal forcing ceases and Ec is conserved in the absence of dissipation. In fact, since 
some numerical diffusion is included, Ec is observed to slowly decay in some cases. This effect is, 
of course, more significant if small scale disturbances are excited as in the case of models with a 
significant solid core. 

A3.2 The canonical Angular momentum 

Similarly, the canonical angular momentum is given by 



Jc = —Xm 



■J V 



(A12) 



where Xm indicates that the imaginary part is to be taken. Here we recall that the state variables 
are complex. This behaves in a manner analogous to the canonical energy, but as applied to the 
total angular momentum content of the planet. 

APPENDIX B: ANALYTICAL EXPRESSIONS FOR THE ENERGY AND ANGULAR 
MOMENTUM TRANSFER DUE TO TIDES ASSOCIATED WITH INERTIAL MODES 

Expressions for the energy, AE^, and angular momentum, AL^, transferred to the planet as a 
result of tides associated with azimuthal mode number m acting during a parabolic encounter 
were derived in PI. For completeness, we briefly review them below. They are: 

y^'^^) k 'I 
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A^2 = j^.^'U-li^ - --DQlil-Jy)}^, (B2) 

where cffc = ak/Vt, with ak being the mode eigenfrequency and the sum is taken over all modes 
with a particular value of m. Contributions from m = and m = 2 have been considered. In 
addition = Vl/Vlp and y = 0(d-fc + m), 

The coefficients C2 = ^, Co = 3/4 and the functions h-miv) are determined by Fourier 
transform of the tidal potential. The latter quantities are described by Press & Teukolsky 1977. 
They are exponentially small for large values of the argument y. 

The overlap integrals Qk measure the coupling of a particular mode with the tidal potential. 
They are given by 

Qk = (^F^'P^l^k) /^k, (B3) 

where Wk is the eigenfunction for mode (k) , N^ = N^/Vl?, P^ is the associated Legendre 
function and F{r) is a correction due to the perturbation of the gravitational potential of the planet. 
In the Cowling approximation adopted here, F = 1. The norm N^ is given by 

Nk = al{Wk\AWu) + {Wk\CWu) (B4) 

where the inner product is defined through 



(l^i|Vr2) = / dz-DjdwWlW2, (B5) 

J y 

with V being the planet volume and the operators A and C are defined below. 



Bl Expressions for the energy and angular momentum transfer due to tides in the 
anelastic impulsive limit 

We here explore the limit in which the characteristic time associated with the tidal encounter ^2^ ^ 
is short compared to the rotation period. In this limit we expect the excitation of inertial modes 
to be impulsive and thus the amount of energy and angular momentum transferred should not 
depend on fine details of the mode spectrum. This is because there is no time for inertial waves 
to propagate, reflect and indicate the presence of any possible standing waves, critical latitude 
phenomena or wave attractors. Later evolution would depend on this but the energy transferred 
would not. This suggests that a relatively simple expression for eg. AE^ should exist which does 
not require detailed knowledge of the mode spectrum. We now show how to obtain this. 

The tidal response to an external potential with fixed frequency a and azimuthal mode number, 
m, satisfies (IP) 
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a^AW-aBW-CW = a^d^{^ext,m-W), (B6) 



where d = 4i7^ — a^, and the operators 

-uo ovo \ aw ) az \ az ) w^ 

B = -^ and C = -4fi2 p . (B8) 

zo ow oz \ az ) 

Adopting a zero density surface, the operators A, B and C are self-adjoint. The eigenfrequencies, 
o = ak and eigenfunctions W = Wk satisfy (IB6I) with \l/ext,m set to zero. 

In the anelastic approximation W is neglected on the right hand side of (IB6I) (IPN) so that the 
tidal response then satisfies 

a^AW - aBW - CW = a^d^^e.t,m. (B9) 

To consider the impulsive limit, we restore the time dependence by replacing a by the operator 
id/dt and note that in the limit of interest where the tidal interaction is fast, a may be considered 
to be large. This suggests we seek a formal solution of (IB9I) as a series in a of the form 

oo 

W = Y,^j^'^'^- (BIO) 

3=0 

Having obtained this we may go on to find the impulsively induced velocity using equation (Hj) 
without viscosity, which leads to 

POO 

Av = - / VWdt. (Bll) 



Here we remark that although we have written infinite limits, we have to suppose that there is 
a formal time scale separation so the limits may be assumed large (and effectively infinite) in 
magnitude when compared to Vt^^ but small when compared to Vt^^ . 

We may then obtain the energy transferred by substituting this velocity into the expression for 
the canonical energy (IA11|) while neglecting the pressure term oc 7~^ which is assumed negligible 
in the anelastic approximation and also the boundary term in that expression. To evaluate (|B11I) , 
recalling the operator interpretation of cr, the assumed time scale separation, and the fact that the 
forcing potential may be assumed to vanish before and after the encounter, one may verify that 
only the term with j = 2 in (|B10I) is needed to evaluate (IBllI) . This is readily found to be given by 

W = J^2 = A-' (4f]2 - CA-i - BA-^BA-i) ■ 4^ext,m (B12) 

Although this expression involves inversion of the elliptic Laplacian like operator A, which de- 
pends on an inner boundary condition that could couple different orders in a way that has not been 
indicated above, this should be in principle straightforward. On following the procedure outlined 
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above to calculate A^E^ (the angular momentum transferred , AL^, is formally higher order in 
the ratio Vt/Vtp so we do not consider it further here) one sees that it scales as Q:^ and that it typi- 
cally involves integral operations on the perturbing potential over the planet volume. Thus a small 
core is expected to have little influence on account of its small volume irrespective of details of 
the spectrum. We finally comment that although the above formal argument required Vt/Vtp to be 
small, on physical grounds we would expect similar conclusions to hold provided there is no time 
for wave propagation with multiple reflections during the tidal encounter. 

APPENDIX C: ANALYTIC SOLUTIONS AND VANISHING OF OVERLAP 
INTEGRALS FOR CORELESS MODELS WITH FIXED QUADRATIC 
GRAVITATIONAL POTENTIAL 

Goodman & Lackner (2009) considered special barotropic models which are in hydrostatic equi- 
librium under a fixed quadratic gravitational potential ^ in the unperturbed state. Thus 

p ar p ar ar 

with \E' = ci + C2r^, with ci and C2 being constant. 

As the Cowling approximation of neglecting the perturbation to the gravitational potential is 
adopted, these models can be treated in the same way as the standard polytropic models we have 
considered which are in hydrostatic equilibrium under a self-consistent gravitational potential in 
the unperturbed state. 

If the equation of state is regarded as being free, a fixed quadratic gravitational potential model 
can be constructed with any density distribution. Thus there will always be one corresponding to 
a polytropic distribution which will have the same anelastic mode spectrum as this depends only 
on the density distribution (see (IB9I) ). However, such models have a special property in common 
with the uniform density model that (l/'Oj)d^/d'cu = {l/r)d'^ /dr = C2 = cOq is constant. A 
consequence of this is that there is an analytic solution that shows that no inertial modes are excited 
in response to forcing by a quadrupole (m = 2) potential as was indeed pointed out by Goodman 
& Lackner (2009) in the case when the fully compressible response governed by equation (IB6I) is 
considered. We here point out that this is also the case when the anelastic approximation is used 
and thus a spurious excitation of inertial modes does not occur in this situation as was suggested 
by Goodman & Lackner (2009). 

To solve the response problem for the quadratic potential models with quadrupole forcing with 
m = 2, we set ^ext,m = Cacc^ where C3 is constant (see equation (IA6l) ) in the governing equation 
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a^AW - aBW - CW = a^d^^^^.t,^ - fW), (C2) 



c 



where we have inserted the constant / which allows us to pass between the fully compressible 
case / = 1 and the anelastic approximation / = 0. A solution can be found by setting W = C/^w'^ 
in (IC2h where C4 is a constant to be determined. If this is done one readily finds that 

2a{(y - 2fi)u;^ + jo'^d) C4 = o'^dc^. (C3) 



When 0" 7^ for any value of / or cr 7^ 2fi for / = 0, this straightforwardly implies that 

C4 = "^^^^ + ^^^ (C4) 

'^ (/a^ + 2n/a - 2u:l) ' "^^^^ 

One may also verify that this solution applies to the cases when cr = and cr = 2fi with / = 
as well, there being zero response to zero forcing in the former case. In the latter case there is 
apparently a unique normal mode with W (xw^ . However, the form of the forcing means that this 
also does not produce a resonant singularity in the response and therefore it is not associated with 
any tidal energy or angular momentum transfer as is also implied by equations (IB II) and (|B2I) . 
In the fully compressible case, / = 1, (|C4I) implies singularities only when o = — f2 ± 



2uq + Q'^ which corresponds to the forward and backward propagating / modes. This means 
that only these modes and no inertial modes would be excited in a tidal encounter. In the anelas- 
tic case, there are no singularities at all meaning no modes are excited. Also, as expected from 
the discussion in IPN, when \a\ and Vl are of the same order and small, the anelastic and fully 
compressible solutions differ by the order of (i^/wo)^- 

It is of interest to connect the above solution to the energy and angular momentum exchanges 
during a tidal encounter given by (IB II) and (IB2I) . From this we expect these exchanges to be zero, 
which in turn implies that for any inertial mode that could contribute, with a^ 7^ 0, or cr| 7^ 40^, 
the overlap integral Qk = 0. It is easy to verify that this is indeed the case by writing down the 
anelastic normal mode equation 

C{Wk) = crlAWk - auBWk - CWk = (C5) 

and considering the inner product 

tu^\C{Wk)) = 0. (C6) 



From this one obtains 
[^\W,) = (C7) 

which confirms that indeed the overlap integrals are all zero when Wq is constant. For the global 
modes that could potentially have significant overlap integrals, adopting appropriate relative er- 
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ror estimates, this result is well represented using spectral approaches of the kind considered in 
Papaloizou & Pringle (1981) and IP as well as the finite difference approach considered in this 
paper. 

We emphasise again that the vanishing of the overlap integrals for the quadratic potential mod- 
els depends on the constancy of cOq. In realistic planetary or stellar models, that are in hydrostatic 
equilibrium under their self consistent gravitational potential, the amount of variation of cuq is mea- 
sured by the central condensation or the ratio of central to mean density which differs significantly 
from unity even for a standard polytrope of index n = 1 as considered in this paper. Accordingly 
there is no reason to suppose that the overlap integrals should be zero in that case. 

This paper has been typeset from a TgX/ ETgX file prepared by the author. 



